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ABSTRACT 

We present a framework for understanding the dynamical and spectral properties of 
X-ray Binaries, where the presence of an organized large scale magnetic field plays 
a major role. Such a field is threading the whole accretion disk with an amplitude 
measured by the disk magnetization (j,(r, t) = B^/(fj, P tot ), where P to t is the total, gas 
and radiation, pressure. Below a transition radius rj, a jet emitting disk (the JED) is 
settled and drives self-collimated non relativistic jets. Beyond rj, no jet is produced 
despite the presence of the magnetic field and a standard accretion disc (the SAD) is 
established. The radial distribution of the disk magnetization [i adjusts itself to any 
change of the disk accretion rate to, thereby modifying the transition radius rj. 

We propose that a SAD-to-JED transition occurs locally at a given radius, in a 
SAD when ji = ji max — 1 while the reverse transition occurs in a JED only when 
fj, = fi m in — 0.1. This bimodal behavior of the accretion disk provides a promising 
way to explain the hysteresis cycles followed by X-ray binaries during outbursts. 

Key words: Black hole physics - Accretion, accretion disks - Magnetohydrodynamics 
(MHD) - ISM: jets and outflows - X-rays: binaries 



1 INTRODUCTION 



In the last ten years, different studies have shown 
that during periods of strong activity XRBs fol- 
low an hysteresis cycle in the hardness-intensity di- 
agram (hereafter HID), the X-ray equivalent of the 
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The archetype of HID has 



been reported in Fig. Q] in the case of GX 339-4 during its 
2002/2003 outburst. The system transits from radio loud 
(jet dominated) hard states (HS) at the right in the figure 
to radio quiet (disk dominated) soft states (SS) at the left. 
Such transitions occur at a luminosity generally a factor 
of a few larger than during the SS-HS transitions back to 
the initial state. This hysteresis cycle clearly indicates that 
the accretion rate cannot be the unique physical parameter 
that controls the spectral evolution of these systems. And 
the fact that at some point XRBs produce collimated jets 
(as probed by flat radio spectra) points towards the other 

physical parameter: the magnetic field. 

In a previous paper (|Ferreira et al.ll2006l hereafter F06), 
we presented a new model to explain the different spectral 



states observed in Black Hole XRB that assumes the pres- 
ence of a large scale vertical magnetic field B z (of the same 
polarity everywhere) anchored in an accretion disk. The cen- 
tral regions have a multi-flow configuration consisting of an 
outer standard accretion disk (SAD) down to a transition 
radius r j and an inner magnetized jet emitting disk (JED) 
below rj. The JED drives a self-collimated non relativistic 
MHD jet surrounding, when conditions for pair creation are 
met, a ultra relativistic pair beam on its axis. We also as- 
sume the presence of a hot thermal corona at the base of 
the j et that is likely part of the JED itself (|Petrucci et alj 
2006, Petrucci et al. 2007 in prep.). This model provides a 
promising framework to explain the canonical spectral states 
of BH XRBs mainly by varying r j and rhQ independently. 
The hard states then correspond to large rj while soft ones 
correspond to small rj. This transition radius is observa- 
tionally identified as the usual "disk inner radius ri n " de- 
termined from spectral fits. Indeed, the JED is transferring 
most of the released accretion power into the jets: the black 
body spectrum is therefore dominated by the outer SAD. 

However this model did not address the cause of the 
transitions between the various spectral states, namely what 
really drives the variations of rj. This directly translates to 
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1 rh = M/M^dd with M the accretion rate and M e dd = ^edd/c 2 
with Lf.dd the Eddington luminosity and c the speed of light. 
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the fundamental question: what are the physical conditions 
to have a JED i.e. to produce a large scale MHD jet from 
an accretion disk? This points out to the important local 
parameter called the disk magnetization 



n(r, i) 



(1) 



where P to t = P gas + Prad includes the plasma and radiation 
pressurefl The disk magnetization is a fundamental param- 
ete r in accretion-ejection phy sics a nd, as sh o wn an alytically 
bv lFerreira fe Pelletierl il 19951) andfFerreiral (|l997l), and nu- 
merically by e.g. Casse fc Keppens 1 20021 ) and Zanni et"ail 
(2007), it has to be close to unity (i.e. the magnetic field 
must be close to equipartition) to produce super- Alfvenic 
jets. We here address the crucial role of the disk magneti- 
zation in the spectral evolution of XRBs. We show that its 
expected variation during the outburst, coupled with that 
of the accretion rate, can explain the hysteresis cycle across 
the HID. 



2 THE DISK MAGNETIZATION 

The reason why the disk magnetization has to be close 
to unity to produce steady-state jets is twofold. On one 
hand, the magnetic field is vertically pinching the accre- 
tion disk so that a (quasi) vertical equilibrium is obtained 
only thanks to the gas and radiation pressure support. As 
a consequence, the field cannot be too strong. But on the 
other hand, the field must be strong enough to accelerate 
efficiently the plasma right at the disk surface so that the 
slow-magnetosonic point is crossed smoothly. These two con- 
straints can only be met with fields close to equipartition. 
Quantitatively, analytical self-similar calculations show that 
fi has to lie between n-min — 0.1 and ji m lLI ~ 1 for a station- 
nary jet from a kepleria n disk to exist (jFerreira fc Pelletierl 
1 19951 . see e.g. Fig. 2 in lFerreirdll997h . The existence of a 
limited range in /i is a key point in our interpretation of the 
hysteresis cycle. 

In our framework, a large scale, organized vertical 
magnetic field B z is threading the whole accretion disk. 
There have been some claims in the literature that such 
a field could not be maintained in a SAD (see e.g. 
iLubow. Papaloizou. fc Pr inglc 1993). In fact, the magnetic 
field distribution is given by the interplay between advec- 
tion and diffusion, namely Ohm's law. It turns out that in a 
SAD in a stationary state, this equation gives B z oc r _7J ™ 
with IZm — He — 1. As a consequence, one can reasonably 
expect n(r) to be a decreasing function of the radius (cf. 
Eq. 4 of F06) . On the other hand, in a JED one has roughly 
B z tx r _5/4 jFerreira fc Pelletier| [l995) and H remains con- 
stant with the radius. In consequence, a reasonable radial 
distribution of /i should look like the one plotted on top of 
Fig. [T] namely constant (but in the range fi m in—Hmax) in 
the JED, located between the disk innermost stable circular 
orbit risco and rj, and decreasing in the SAD beyond rj. 

Now, the evolution in time of the disk magnetization 
fi(r,t) is a far more complex issue. Since Ptot oc m, one gets 



2 Note that fi differs from the inverse of the usual plasma /3 pa- 
rameter which compares the magnetic pressure to the gas pressure 
only. 



fi cc m~ 1 B z . It means that its temporal evolution is gov- 
erned by the evolution of both B z and m(t). Although the 
origin of the variations of the disk accretion rate rh (t) dur- 
ing a n outburst is relatively well understood (e.g. lLasotal 
l200ll ). the evolution of B z (t) is more speculative. One of 
the reasons is that B z could itself be a function of the disk 
accretion rate. Different scenarios h ave been proposed (e.g. 
iTagger et aill2004l ; iKing et al.ll2004 ) but the problem is very 
tricky and no detailed calculation has emerged so far. In 
this context, and given the large uncertainties on the phys- 
ical processes governing the evolution of the magnetic field 
with m, we make the zeroth order approximation that the 
variation of the accretion rate is the dominant cause of the 
variation of [i. In other words we assume here that fl(r,t) 
anti-correlates with the accretion rate. 

One last issue is what determines the transition radius 
r j in a disk that is threaded everywhere by a large scale B z 
field? In a SAD, the origin of the turbulence is commonly 
believed to arise from the magneto-rotational instabilit y 
or MRI (jBalbus fc Hawlevlll99ll ; Uesur fc Longarettill2007f l. 
However the MRI requires the presence of a weak magnetic 
field and is quenched when the field is at equipartition. Thus, 
while a SAD works well at /i <C 1, there is no known solution 
when fi > 1. On the other hand, as said before, a steady- 
state JED can establish only around a rather limited range 
in disk magnetization, namely fimin < /i < fimax- 

These facts lead us to make the following conjecture 
that, at a given radius, 

- a SAD-to-JED transition occurs when fj, = fi ma x — 1, 

- a JED-to-SAD transition occurs when /j, — fi m in —0.1. 
As we will see next Section, this conjecture provides a 

plausible scenario for the hysteresis cycle observed in BH 
XRBs. 



3 TEMPORAL EVOLUTION OF BH XRBS 

It is commonly believed that BH XRBs generally evolve 
along an hysteresis cycle during their outbursts, following 
a counter clockwise sequence A-B-C-D like the one reported 
in Fig. [T] before turning back to their initial state at the 
end of the outburst. The HS and SS correspond to the A-B 
and C-D parts of this diagram respectively. The horizon- 
tal branches, B-C and D-A, correspond to the intermediate 
states (hereafter IS) whos e spectral and timing properties 
are relatively comp l ex (e .g. iBelloni et al.|[2005l ; iFender et alj 
|2004 ICorbel et alJliooi 1 ). This is especially the case dur- 
ing the high IS (hereafter HIS, iHoman fc Bellonil 120051 ; 
IBelloni et all [2005) that corresponds to the top horizontal 
branch where the system luminosity is high. During its evo- 
lution along this horizontal branch the system may enter in 
a flaring state associated with radio an d X-ray flares and /or 
superluminal sporadic ejections (e.g. ISobczak et alj |2000| ; 
lHannikainen et all 1200 ll ; Vender et all \2004 ). This~spectral 
evolution goes also with a quenching of the radio emission, 
signature of the disappearance of the jet component. Dur- 
ing the bottom horizontal branch (the low IS, hereafter LIS) 
and the turning back of the system to its initial state, the 
radio emission is detected again, meaning that the ejection 
structure has been rebuilt. 

We detail below our interpretation of this hysteresis cy- 
cle as a consequence of our conjecture on the disk magneti- 
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Figure 1. Left: Light curve and Hardness-Intensity Diagram followed by GX 339-4 during its 2002-2003 outburst (from iBelloni et all 
2005). Right: Sketch of the expected evolution of the radial distribution of the disk magnetization £t(r) during an outburst with hysteresis. 
In each plot, the solid line represents the initial (for each part of the HID) radial distribution of the disk magnetization, while the dashed 
one shows its evolution. See Sect. [3] for more details. 



zation fj,(r,t). We give also some estimates of the expected 
time scales involved at the different phases of the cycle. In 
the accretion disk, the shortest dynamical time scale dealt 
with is the keplerian orbit time, namely td{v) = 2-K/Q.k = 
0.1 (m/10)(r/50rg) 3//2 sec, where m is the black hole mass 
in solar mass units and r g = GM/c 2 the gravitational ra- 
dius. In the following we take m — 10 and consequently 
r g = 1.5 x 10 6 cm. 

Let us start at a Low/hard State located at the bot- 
tom of the HID right branch (point A in Fig. [T|. This 
state is usually seen as a SAD being interrupted at typi- 
cally ri n ~ 100r 9 . Within our framework, this is perfectly 
consistent with a JED settled from risco up to rj = ri n 
and driving self-collimated non relativistic jets (seen as per- 
sistent radio features). Within the JED \x is roughly con- 
stant and we assume that at the beginning of the outburst 
H = Umax — 1. Let us now consider an increase in the disk 
accretion rate rh (an outburst) triggered by, for inst ance, a 
disk instability at some outer radius r out 3> Tj (e.g. lLasotal 
|2001| ). A front defined by a rise in P to t (hence a decrease 
in /i) propagates then radially inwards eventually reaching 
r.j. As for the SAD, the JED adapts itself locally on a dy- 
namical time scale (a few td) by lowering fi. This adap- 
tation will affect the whole JED extension on a timescale 



roughly equal to T acCt jED, the accretion time scale in the 
JED. Si nce the accre tion radial velocity w' ED ~ Q.k h in 
a JED (|Ferreiralll997l ). this time is even shorter than that 
of the front propagation in the outer SAD. For instance, 
one gets r acc , jED(W0r g ) ~ T D (100r a )r/h ~ 30 sec (usin g 
a disk aspect ratio h/r = 0.01. IShakura fc Sunvaevlll97ot ). 
This timescale is much shorter than the timescale of the as- 
cendent phase of an outburst, which is typically of the order 
of several days between states A and B. Thus, during the 
ascendant motion of the XRB along the right branch from A 
to B, /1 should only decrease within the whole JED without 
any change in rj: no spectral state transition is expected. 

If the rise in m is such that /i(rj) decreases below the 
critical value fimin, then the XRB arrives at point B. Any 
subsequent increase of the accretion rate will then produce 
an outside- in JED-to-SAD transition starting from the outer 
part of the JED. The transition radius rj decreases and both 
the MHD jet and the corona extension shrink. The XRB 
enters the high intermediate state, between B and C. Since 
r j is decreasing, the timescale for the system evolution along 
the top branch become s more and more ra pid, in agreement 
with observations (e.g. IBelloni et al]|2005l ). 

Formally the end of the top branch should be reached 
around point C when rj = risco i.e. when the JED has 
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completely disappeared and the whole disk adopts a radial 
structure akin to the standard disk model. However, we ex- 
pect this situation to be quite variable, probably even before 
rj — risco- Indeed, an important large scale (open) mag- 
netic field is present, providing a magnetic torque helping 
the MRI to drive accretion. But it is not clear whether jets 
can be formed or not. In fact, the only calculations that ac- 
tually proved that steady ejection can take place from a disk 
were done with a self-similar ansatz. Implicitly, this assumes 
that the radial extension of the jet emitting region is large. 
This is consistent with the physical interpretation that the 
disk material must be first radially accreting fast before be- 
ing lifted up by the plasma pressure gradient and magneto- 
centrifugally accelerated at the disk surface by the Lorentz 
force. What will happen when rj > risco only? Our guess 
is that the MHD jet structure will probably be oscillating 
between life and death. Accordingly, persistent radio jets 
should thus disappear even before point C is reached. Be- 
sides, in our framework, we also expect luminous systems to 
enter a flaring state as soon as the conditions for a strong 
pair production in the MHD funnel are fulfillecH resulting in 
sporadic ultra-relativistic blob ejections. This could explain 
why this HID r egion seems to harbor complex variab ility 
phenomena (e.g. iBelloni et al.ll2005l : iNespoli et al.ll2003h . 

In C, the system enters the so-called Soft State. The 
disappearance of the MHD jet would explain the quench- 
ing of the radio emission generally observ ed at this stage 
(e.g. lCorbel et al.ll2003l : IFender et al.ll2004l ). Then whatever 
the accretion rate does, and following our conjecture, the 
JED cannot reappear along the left vertical branch unless 
h{tisco) reaches fi m ax ~ 1. The descent from C to D must 
then correspond to a decrease in intensity i.e. mainly to a 
decrease of the accretion rate itself, with no spectral state 
transition. This is the beginning of the fading phase of the 
outburst. In our simple picture, the decrease of the disk ac- 
cretion rate produces an outside-in decrease in P to t leading 
to an increase in fi. Point D then corresponds to the situa- 
tion when n{risco) = /J-max- 

This initiates the inside-out rebuilding of the accretion- 
ejection structure along the low horizontal branch. Note that 
the accretion rate must still be decreasing. But it does not 
require a strong decrease in rh to build up again a JED 
on a significant extension (hence jets observable in radio 
and a change in the hardness ratio). Indeed, as said before, 
the steady-state radial distribution of the vertical magnetic 
field is steeper in a SAD (B z oc r _7Jm , with lZ m ~ 1) than 



in a JED (B z oc 



-5/4\ 



As a consequence, \i may become 



greater than unity in the innermost disk regions with only 
a slight decrease in rh. There is no known steady state disk 
solution at these large values of magnetization. We expect 
such a strong field to diffuse away outwardly in the outer 
SAD so that a maximum value of /x = 1 can be maintained 
in the inner JED-like regions. Providing a time scale is out 
of the scope of the present paper as it would require to 
take into account the whole process of the magnetic field 
redistribution. However, this time is probably shorter than 
(in the JED) or comparable to (in the SAD) the accretion 
time scale. We thus expect a rather fast inside-out rebuilding 



3 Thi s would correspond to the jet line defined bv lFender et al.l 
J2004h , see F06 for more details. 



of the JED with reappearance of the self-collimated electron- 
proton jet and the hot corona. This marks also a spectral 
state transition as the SAD emission is receding when going 
back to point A, at the end of the outburst. 



4 DISCUSSION AND CONCLUSION 

We present a scenario where the spectral states of X-ray 
binaries depend on the transition radius r.j between an inner 
accretion-ejection structure and an outer standard accretion 
disk (F06). We propose in this paper that the variation of 
rj is mostly controlled by the evolution of the accretion rate 
rh and the disk magnetization /i. 

We made the conjecture that a SAD-to-JED transi- 
tion will locally occur when fi reaches a maximum value 
Umax — 1, whereas a JED-to-SAD transition requires a 
much lower fi : analytical estimates show that it could hap- 
pen around fi m in — 0.1. It was shown that such a simple 
framework is enough to explain the dynamical and spec- 
tral behaviors of XRBs, providing moreover a dynamical 
basis for the hysteresis cycles observed in the HIDs. The 
exact values taken for the limits [imin and fi m ax may not 
be strictly equal to those cho sen here as th ey were derived 
from self-similar calculations <|Ferreiralll997h . However, this 
is not crucial as the main element is the existence of such 
an interval. It is noteworthy though that these values are 
compatible with our current understanding of disk physics: 
there are neither SAD solutions with large scale fields be- 
yond equipartition, nor JED solutions at much lower fields. 
However, a precise examination of the dynamical transitions 
between a JED to a SAD solution and vice versa remains to 
be done in order to assess the conjecture made here. 

Interestingly this scenario gives a physical constraint 
on the existence of a hysteresis cycle in a given source dur- 
ing an outburst. Indeed, the magnetization has to decrease 
necessarily below /J, m in in the whole accretion disk during 
the flare. This translates into the condition that, at a given 
time during the burst, n(risco) has to become smaller than 
l-imin ■ If it is not the case the source should oscillate between 
hard and soft states with an inner JED never disappearing. 
Cyg X-l could be in this situation. It is known to transit 
from hard to soft state and reversely at the same accre- 
tion rate rh ~ 0.01 (e.g. IZdziarski et alJliool I Smith et al.l 
2002). It is persistent in X but also in radi o, even if the ra- 
dio e mission is weaker in its softest states ( Brocksopp et al.l 
1 19991 and references therein, iPandev et alj 20061 ). This sug- 
gests that ejection is always present. On the other hand the 
strong disk signature in the soft state implies that rj is 
close to risco- Thus for some reasons, in this object the 
accretion-ejection structure keeps stable even for low value 
of rj. That could be linked to the fact that it is a high mass 
X-ray binary system as alr e ady p roposed by ISmith et al.l 
l|2002l ); lMaccarone fc Coppl l|2005T) . 

GRS 1915+105 has also a peculiar behavior that does 
not match the one described in Sect. [3] This source is gen- 
erally believed to accrete close to its Eddington limit and to 
stay preferentially at the top of the HID. Its "hard" states 
would correspond to the right part of the top branch while its 
soft states are generally observed i n the left part of th e top 
branch, after the so-called jet line (iFender et alj |2004). The 
flare states are observed when the source oscillates around 
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this line. In our view, in this region of the HID the JED com- 
petes with the SAD with rj ~ risco- The magnetic field 
keeps trying to re-adjust itself in the inner disk region to 
agree with one of the two disk solutions. This might produce 
the complex va riability pattern exhibited by this source dur - 
ing outbursts l|Belloni et al.ll2000l ; iFeiider fc Belloni l2004h . 
Moreover, close to the ISCO the time scales controlling the 
system are of the order of fractions of second, in agreement 
with the rapid variability that characterizes this peculiar 

object. 

Maccarone & Coppi (2003) and Zdziar ski fe Gierl inski 

(2004) provided several other suggestions for the hysteresis 
behaviour of XRBs but without going much into de tails (see 
critical discussion in iMever-Hofmeister et al"1l2005l ). In fact, 
it seems that t he most detailed model so far is the disk evap- 
oration model ( Mover- Hofmeister et al. 2005 and references 
therein). In this model the efficiency of the SAD evaporation 
critically depends on the external photon field. In the Soft 
State for instance, there is mainly a Compton cooling and 
not much evaporation. The spectral state transitions would 
then be triggered at different accretion rates simply because 
of history, hence of a different Comptonizing spectrum (see 
their Fig. 2). Although the physics of jet production is not 
included in this model, taking into account radiative effects 
in order to explain the HID hysteresis sounds quite reason- 
able and could, in principle, be integrated in our framework. 
In our case however the hysteresis has a dynamical explana- 
tion based on the disk magnetization and a two-temperature 
disk never settles in. Since ejection processes are known to 
be important in the XRB global energy budget, our suspi- 
cion is therefore that evaporation is not the main driver. 

As a final note, we stress that if one is to explain self- 
collimated jets in XRBs, then large scale magnetic fields 
in the inner regions of accretion disks must be taken into 
account. Their presence is both compatible with standard 
accretion disk physics and provides a simple and promis- 
ing framework to explain hysteresis cycles. In our view, the 
magnetic flux available in disks is therefore a fundamental 
and unavoidable ingredient for modeling XRBs. Also, there 
is no reason why it should not vary from one system to an- 
other. Since changing the amount of magnetic flux changes 
the transition radius rj, the characteristic value of m (hence 
luminosity) associated with each spectral state transition is 
also modified. This picture becomes even more complex if 
accreting material is advecting magnetic flux from the sec- 
ondary. Indeed, taking into account the advection of a large 
scale magnetic fields within the disk introduces a whole new 
set of variable phenomena. 
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